Coherent acceleration of Bose-Einstein condensates 
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We present a theoretical analysis of the coherent accel- 
eration of atomic Bose-Einstein condensates. A first scheme 
relies on the 'conveyor belt' provided by a frequency-chirped 
optical lattice. For potentials shallow enough that the conden- 
sate is not fragmented, the acceleration can be interpreted in 
terms of sequential Bragg scattering, with the atomic sample 
undergoing transitions to a succession of discrete momentum 
sidemodes. The narrow momentum width of these sidemodes 
offers the possibility to accelerate an ultracold atomic sample 
such as e.g. a Bose-Einstein condensate without change in 
its momentum distribution. This is in contrast to classical 
point particles, for which this kind of acceleration leads to a 
substantial heating of the sample. A second scheme is based 
on the idea of a synchronous particle accelerator consisting of 
a spatial array of quadrupole traps. Pulsing the trapping po- 
tential creates a traveling trap that confines and accelerates 
the atomic system. We study this process using the concept 
of phase stability. 

PACS numbers: 03.75.Fi, 32.80.-t, 42.50.Vk 



I. INTRODUCTION 

A number of applications of atom optics, from inertial 
sensors to nanofabrication, atomic holography (^-^] 
and certain schemes for quantum computation , can 
benefit substantially from manipulating quantum degen- 
erate atomic beams with finely controlled potentials. In 
particular, the atomic de Broglie wavelength can be easily 
modified by accelerating or decelerating an atomic beam 
to the appropriate velocity, thereby making it a tunable 
source. There is also an increased need to control ultra- 
cold samples of atoms in confined geometries, as e.g. in 
atomic waveguides for potential integrated atom optics 
applications. First transport experiments in waveguides 
were recently carried out with non-condensed atomic 
samples, both in stationary and in time-dependent po- 
tentials 

The simplest and most important coherence property 
of an atomic beam is its monochromaticity. With this in 
mind, our goal in this paper is to determine the extent to 
which it is possible to significantly accelerate an initially 



quasi-monochromatic atomic sample without losing this 
property. For this purpose it is sufficient to describe the 
atomic system in the Hartree mean- field limit, remember- 
ing that this "classical-like" description makes implicit 
assumptions on the quantum statistical properties of the 
Schrodinger field, essentially assuming that it is in a co- 
herent state. Other, more subtle aspects of the quantum 
coherence of accelerated matter-wave fields JTl|~p^ will 
be addressed in future work. 

We consider two different acceleration schemes: in the 
first one, an accelerated trap is provided by a frequency- 
chirped optical lattice, while in the second one we con- 
sider a synchronous accelerator consisting of a spatial 
array of quadrupole traps, which is similar to a scheme 
that was successful in manipulating the velocity of neu- 
tral molecules and atoms |l^JT^. We restrict our discus- 
sion to one dimension for simplicity. 

Ultracold atoms in both static and time-dependent op- 
tical lattices have been investigated in several contexts in 
the recent past. At the single-atom level, they were ex- 
ploited extensively in theoretical and experimental work 
aiming at demonstrating effects such as Bloch oscillations 
]r7[-p^ , Landau-Zener tu nnehng l2 C|j2l}| , the appearance 



of Wannier-Stark ladders 




quantum chaos |2J 
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and the dynamics of mesoscopic quantum superpositions 
p6t . The extension of this work to Bose-Einstein conden- 
sates |2^,^ includes the first demonstration of a mode- 
locked atom laser a device that can also be inter- 
preted in terms of Landau-Zener tunneling. Bragg scat- 
tering of Bose-Einstein condensates off a periodic optical 
potential was also used in experiments to characterize 
their coherence and phase properties [^^. It was ex- 
ploited to generate multiple momentum sidemodes of the 
condensate used in matter-wave four-wave mixing exper- 
iments [^-|4j and to realize atomic beam splitters used 
in coherent matter- wave amplifiers | p5| , p6| . 

Of particular relevance in the present context is Ref. 
p9[ , which explicitly considers the acceleration of atoms 
in a moving periodic potential and interprets it in terms 
of momentum transfer via multiple adiabatic rapid pas- 
sage. This reference also reports experiments demon- 
strating the acceleration of atomic molasses with large 
initial momentum spreads as well as of ensembles with 
subrecoil momentum distributions. The final momen- 
tum spread is limited by the the lattice recoil energy for 
the broad samples and by the initial distribution for the 
subrecoil atoms, see Figs. 11 and 13 of Ref. |l|. We 
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will show that we can suppress heating even in the case 
the atomic temperature is substantially lower than the 
recoil temperature associated with the lattice, which in 
particular holds for a Bose-Einstein condensate and its 
quasi-monochromatic momentum distribution, n 

The paper is organized as follows: Section |^ investi- 
gates the acceleration of an ultracold atomic ensemble in 
a frequency-chirped optical lattice. We set the stage by 
first considering an ensemble of classical point particles 
spread over a large number of lattice periods. We com- 
pare its dynamics to that of a Bose-Einstein condensate 
of the same spatial extent, described in the Hartree mean- 
field regime. We demonstrate that while the accelerated 
classical particles undergo considerable heating, which is 
however substantially eliminated by the quantum inter- 
ferences associated with the extended atomic wave func- 
tions of ultracold samples. The quantum dynamics are 
interpreted as continuous Bragg scattering in a lattice 
with time-dependent detuning. We briefly review the 
formalism of Bragg scattering within the framework of 
a partitioned momentum space et al. | |3^ , which allows 
one to discuss the acceleration in terms of a coupled- 
mode approach that leads to a determination of opti- 
mum acceleration parameters. We conclude this section 
by some comments on the role of the mean-field interac- 
tion. Section III then turns to the synchronous accelera- 
tion scheme. Again, we first consider the acceleration of 
non-interacting classical point particles, this time using 
the concept of phase stability. We then study numeri- 
cally the acceleration of a condensate. Finally, section 
[V is a conclusion and outlook. 



5{t) = uj_ {t) - uj^ {t) = c [k^ {t) - (t)] 



(3) 



II. ACCELERATED OPTICAL LATTICE 

A. Classical acceleration 

To set the stage for our discussion, we first consider the 
acceleration of an ensemble of non-interacting, point-like 
classical particles in a frequency-chirped optical lattice 
formed by two counterpropagating laser beams. The re- 
sulting time-dependent optical potential is 



where 



V{z, t) = Vb cos[K[t)z - 5{t)i\, 



K{t) = k^{t) + k^{t) 



(1) 



(2) 



is the sum of the instantaneous wave numbers of the left- 
and right-propagating laser fields forming the lattice, and 



Shortly after submission of this work, this prediction was 
independently confirmed in experiments on the acceleration 
of a condensate by a frequency-chirped optical lattice ^] . 



is the instantaneous frequency detuning between them. 
The instantaneous phase velocity of the lattice fringes is 
therefore 



(4) 



which reduces to the usual phase velocity 5/K in the case 
of constant detuning. Assuming that i^iat remains much 
less than the velocity of light c at all times we have 



k^{t) - k^{t) < k^it) + k^{t), 



(5) 



and it is an excellent approximation to set k^{t) + 
k^{t) w 2kL, where k^ is a nominal laser wave num- 
ber independent of time. The potential (|l|) reduces then 
to 



V{z, t) « Vo cos[2kLZ - d{t)t] 



and 



Wlat (t) 



1 



dt 



(6) 



(7) 



where the time derivative is the instantaneous frequency 
of the field. While we have shown numerically that this is 
not necessarily an optimum choice, we restrict our anal- 
ysis in this paper to linear accelerations produced by a 
time-dependent detuning of the form 



Sit) = 7jt, 

producing a lattice group velocity 



(8) 



(9) 



and a constant acceleration aiat = v/^l- This simple case 
is sufficient to identify and discuss the major physical 
mechanisms at play in the acceleration of the atoms. 

To determine the dynamics of an ensemble of classi- 
cal point particles on this lattice, we have numerically 
solved Newton's equations of motion for a large number 
of atoms initially at rest at random points in an interval 
large compared to the lattice period. The size of this in- 
terval and the probability of finding an atom at a given 
point were chosen to mimic the shape of a Bose-Einstein 
condensate density profile, to which we turn shortly. 

The solid lines in Figs. and I show the evolution 
of the mean position {z{t)) and mean velocity {v{t)) of 
the classical atomic ensemble over the acceleration time 
for two values of aiat . We observe that the mean atomic 
velocity is always somewhat less than the instantaneous 
lattice velocity wiat, the dashed line in Fig. ||. The cause 
of this difference is revealed in Fig. |[ which shows the 
momentum distribution of the classical ensemble n{k) at 
a fixed time t and for the same two values of aiat, where 
we expressed the momentum p in terms of the wavenum- 
ber k — p/Ti. In addition to peaks at positive momenta 
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indicative of accelerated atoms that contribute to an in- 
crease in the mean velocity of the sample, a significant 
group of particles acquire negative momenta, i.e. they 
are accelerated in the direction opposite to the lattice 
motion. These are atoms that spill into a well to their 
left in the moving potential. A similar effect will also 
be encountered in the synchronous particle accelerator 
of section III . 

Another important feature of Fig. ^ is that the mo- 
mentum distributions of both the accelerated and decel- 
erated groups of atoms are rather wide. This shows that 
(even ignoring the decelerated atoms), the lattice acceler- 
ator produces a considerable heating of the atomic sam- 
ple. This can potentially be a very serious problem for 
applications requiring a high degree of spatial coherence 
of the atomic beam. Surprisingly perhaps, we will see 
that the situation can be significantly improved in the 
quantum regime, a result of quantum interferences. 



B. Condensate acceleration 

We now turn to the acceleration of a Bose-Einstein 
condensate at temperature T = 0. For the present study, 
it is sufficient to describe its dynamics in a Hartree mean- 
field approximation via the Gross-Pitaevskii equation for 
the normalized condensate wave function ^'(z, t), 



-I 

- — -^ + Vocos[2kLZ-Sit)t] 
NUom''-i'{z,t), (10) 



where N is the number of atoms in the condensate. For 
the reduction to one dimension we assume the conden- 
sate to be in its transverse ground state Ug{rj_). This is 
accounted for by the effective nonlinear parameter Uq, 



(11) 



where Us is the s-wave scattering length. 

In the presence of a periodic potential, one needs to 
beware of the possible fragmentation of the condensate 
via a Mott insulator transition. This problem has been 
investigated in Ref. p9t , which discusses the transition 
to this phase from the superfluid phase as a function of 
the ratio k between the depth of the optical lattice and 
the inter- well tunneling matrix element. The first Mott 
transition, corresponding to an occupation of one atom 
per lattice site, occurs for n « 5.8z, z being the number 
of nearest neighbors corresponding typically to lat- 
tice potential depths of the order of 2.5-5 lattice recoil 
energies. Here we consider only situations where the lat- 
tice potential is shallow enough that such fragmentation 
does not occur. 

We solve Eq. (|o|) numerically using a split-operator 
technique for a condensate of the same initial density 
profile as the classical ensemble considered earlier. The 



resulting mean position {z{t)) and mean velocity {v{t)) 
of the condensate are shown by the dotted lines in Figs. 
|l| and ^. We observe that whereas the displacement dy- 
namics is very similar to the classical case, the mean 
velocity exhibits two differences. First, the condensate 
shows pronounced time-dependent oscillations, the well- 
known Bloch oscillations. Second, the higher of the two 
lattice accelerations causes (v) to saturate to a constant 
value, an effect due to increased Landau-Zener tunnel- 
ing.0 

The difference between the classical and quantum sit- 
uations is evidenced even more strikingly in Fig. IJ. In- 
stead of being composed of two broad continua as in Fig. 
^, the condensate momentum distribution ^!)(fc), which is 
the Fourier transform of the spatial wavefunction ^'(z), 
consists of a series of very narrow peaks located at integer 
multiples of 2fci.0 Physically, this is due to the fact that 
while the classical particles probe the local value of the 
lattice potential, the ultracold Bose-Einstein condensate 
probes a region of the lattice comprising a large number 
of maxima and minima, and hence its full periodicity. 
As such, the individual classical particles are channeled 
in one or the other potential well of the time-dependent 
potential, while the ultracold atoms are diffracted by it. 
Their dynamics is governed by the quantum interferences 
that give rise to Bragg scattering, with the resulting nar- 
row peaks of Fig. ^ This is a central result of this paper: 
it shows that after acceleration, the ultracold atoms can 
still be largely monochromatic. The next section takes 
advantage of the physical difference in physics between 
the acceleration of quantum and semiclassical samples to 
investigate ways to achieve this goal. 



C. Sequential Bragg resonances 

With atom interferometric applications in mind, we 
call an acceleration scheme ideal if it leaves all atomic 
population in one single momentum mode whose value is 
determined by the velocity wiat of the optical lattice. We 
can therefore define a figure of merit of the accelerator as 
one minus the fraction of atoms in other momentum side- 
modes. Fig. ^ suggests that we have to restrict ourselves 
to certain acceleration rates and times in order not to lose 
too many atoms. The closest we were able to approach 



^Superfluid effects are not an issue in our system: The crit- 
ical velocity for typical Bose-Einstein condensates was re- 
cently determined to lie in the mm/s regime pc|-^, much 
less than the velocities to which the condensate is acceler- 
ated. For the acceleration rates under consideration here, the 
instantaneous lattice velocity exceeds the critical velocity al- 
ready after a small fraction of the total acceleration time. 

■^The experiment of Ref. shows clear evidence of such 
peaks. 
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"perfect" acceleration with our simple linear acceleration 
scheme is the example of Fig. §(a). 

In the following we consider a simple coupled-mode 
description of the acceleration process to identify the im- 
portant parameters in its optimization. We neglect the 
effects of the mean-field energy for now and extend the 
case of a static lattice as shown in Ref. to accelerated 
ones. Starting from the coUisionless form of the Gross- 
Pitaevskii equation ( |lO| ) and introducing the momentum 
space condensate wave function that we already used in 
Fig. I 



and the coupling matrix V(t) is given by 



,^(fc,t) 



/27r 



dz '>p{z, t)e 



-ikz 



(12) 



one finds readily that its evolution is governed by the 
coupled difference-differential equations 



d 



2m 
2 



-m 



^From these equations, it is clear that states with mo- 
mentum k only couple to neighboring states with k = 
±2kL- Thus we can partition the momentum space into 
intervals of width 2fci and substitute for (j){k,t) a set of 
wave functions restricted to these intervals, 

(t)^{k) = (t){k) for {n-l)kL < k < {n + l)kL. (14) 

The partitioned wave functions (pn are then expressed 
in terms of the shifted momentum q — k — 2nfci and 
multiplied with an interval-dependent phase factor, 



0„(<z) = 0„(fc)e-""'*'. 



(15) 



Finally, we use the numerically established property that 
each of these partial wave functions suffers a momentum 
spread small compared to the mode spacing , see Fig ^, 
to approximate them by plane waves at g = 0, 



,(q = 0,t) = c„(t). 



(16) 



In units of the lattice recoil frequency t^rcc — 2Tik\/m, 
we introduce the dimensionless variables 



V = V/^rcc, Vo = Vo/huj^cc, (17) 



in terms of which the equation of motion Eq. ( [T^ ) can 
be rewritten in matrix form as 



d 



(18) 



where we defined the vector of coefficients — which es- 
sentially gives the probability amplitudes for the various 
momentum sidemodes — as 



c(t) = (•••,c_i,co,ci,---) 



(19) 



V(r) 



■■ W-i(t) 

^ 





V 




Vo 
2 

wo(t) 

% 
2 







Vo 
2 

uji{t) 



^ 



(20) 



The diagonal elements of V{t) are responsible for the 
time-dependence of the coupling matrix and are given 

by 



(21) 



Finally, we can define the detuning A„ between adjacent 
modes as 



A„(r) = w„(r) - w„_i(t) = 2n - 2f]T ~ 1. 



(13) 



(22) 



Equation ( |18| ) describes Bragg scattering of atoms 
off the periodic optical lattice. Bragg resonances oc- 
cur whenever one of the detunings A„ (r) becomes equal 
to zero. From Eq. (|2^), we observe that these detun- 
ings depend both on the scattering order n and on time, 
this latter dependence resulting from the acceleration of 
the optical potential. As a result of the linear acceler- 
ation, neighboring pairs of modes are therefore succes- 
sively moved in and out of resonance, so that in contrast 
to the case of classical particles, the physical process un- 
derlying the atomic acceleration are successively tuned 
and detuned Bragg resonances.^ 

This sequence of resonances is illustrated in Fig. ||, 
which shows the evolution of the population dynamics of 
a few momentum sidemodes of the condensate. The solid 
lines give the results of a truncated coupled-modes analy- 
sis, while the dotted curves show the results of the direct 
solution of the Gross-Pitaevskii equation, in which case 
the various sidemode populations are calculated from 



(2Tl+l)fc£ 



dk |0(fc,r)|' 



(23) 



{2n-l)kL 



The two approaches are in excellent agreement, despite 
the fact that the coupled-mode analysis included only 
eight modes in the present example. The upper graph 
clearly illustrates the sequential population transfer to- 
wards sidemodes of higher momentum. The lower graph 
focuses on the first side mode of negative momentum. Af- 
ter initial oscillations, indicative of its coupling to other 



*Ref. refers to this process as " multiple adiabatic rapid 
passage," an equally appropriate terminology. 
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modes, its population stabilizes to a final value resulting 
from the fact that it is then far-off resonance from any 
other mode. 



D. Optimal acceleration 

Having identified sequential Bragg resonances as the 
physical mechanism of acceleration of ultracold atoms, 
we can now attempt to optimize the external parameters 
to achieve maximum velocity and figure of merit of the 
accelerator. Ignoring for now all but the two adjacent 
modes c„ and c„_i, we observe that they couple with a 
time-dependent Rabi frequency f2„ given by 

^^n(r)-i\/fo' + A2(r). (24) 

We can gain some insight into the time scale of this cou- 
pling mechanism by introducing an averaged Rabi fre- 
quency f2 over the interval th = 1/fj, the time it takes 
the system to move from one Bragg resonance to the 
next, see Eq. (|2^), 

Tn+TR 

n=— [ dT^nir), (25) 
tr J 

where Tn — {n~ l)/fi. Note that Cl is independent of the 
index n. Hence, one can expect that a close to optimal 
mode-to-mode coupling should correspond to 

Qtr w tt, (26) 

since in this case the system can complete half of an 
averaged Rabi cycle in the time it takes to move from 
one Bragg resonance to the next. The complete popu- 
lation of one mode is then approximately transferred to 
the next mode, thereby increasing the momentum of the 
condensate by 2hk]^. If the lattice acceleration is too fast 
for this Rabi transfer to fully occur, lower momentum 
modes remain significantly populated. This is the case 
in the example of Fig. ^(b) . Conversely, atoms acceler- 
ated too slowly undergo more then one half Rabi cycle 
between two modes. This results in oscillations in the 
expectation values (v) and (z). 

We have mentioned in the context of Fig. ||(b) that 
once a momentum sidemode is shifted out of resonance, 
its population remains constant. This accounts for the 
saturation in (v) shown in Fig. || In contrast to the clas- 
sical case where all atoms that are captured in a poten- 
tial well remain captured, here, in the quantum regime, 
atoms initially accelerated may gradually be moved out 
of resonance, after which they retain a constant velocity. 

Fig. 1^ summarizes the results of a numerical optimiza- 
tion of the lattice acceleration rate. It confirms that if 
the acceleration rate is too large one can only transfer a 
small fraction of the population to the next mode, while 



if it is too small we couple to modes with negative mo- 
mentum, as manifested in the small dip in all curves near 
the origin. The plateau-like feature defines the regime of 
efficient coupling, where around 85% of the population is 
transferred. 

We mention that it is also possible to accelerate the 
condensate by applying a series of Bragg pulses, each of 
them resonant with the next higher pair of momentum 
modes. For the specific parameters used in our figures, 
one could then reach accelerations approximately twice 
as high as in the case of continuous detuning, but further 
improvement would require stronger Bragg pulses that 
create additional couplings to different modes |3^]. We 
do not go into further details about this approach here. 
Rather, future work will implement learning algorithms 
to determine ideal chirping and/or pulse sequences to 
achieve the coherent and lossless acceleration of ultracold 
atomic samples. 

We conclude this section by remarking that so far, we 
have ignored the effects of the mean-field energy on the 
Bragg acceleration of the atomic sample. But it should be 
expected that the associated nonlinear phase shifts can 
effect it significantly. This is illustrated in Fig. |7[ where 
we plot the temporal evolution of the zeroth momentum 
sidemode for different mean-field energies (atomic num- 
bers). For large mean-field energies we observe that the 
local evolution can deviate significantly from the low- 
density limit, although the general behavior is similar. 
For the Sodium parameter used in our simulations, we 
find that in condensates of up to 10^ atoms the nonlin- 
ear effects are small enough that the system obeys the 
predictions of the simplified linear model. 

III. SYNCHRONOUS PARTICLE 
ACCELERATION 

We now turn to our second accelerator scheme, which 
uses a time-varying conservative potential well, formed, 
e.g., by external magnetic fields. The idea is to initially 
place the condensate near the top of a potential well. As 
it moves down the well, it gains kinetic energy at the 
expense of its potential energy. If the external field is 
turned off precisely when the atoms are near the poten- 
tial minimum, then the condensate will wind up with an 
increased kinetic energy. The process can be repeated by 
letting the condensate pass through a sequence of pulsed 
potential wells appropriately displaced in space and time. 
It can thus gain a considerable amount of kinetic energy. 

We first discuss how such a method can be used to 
accelerate an ensemble of classical particles. Fig. || illus- 
trates the potential energy, V{z), as a function of particle 
position. In analogy to the concepts of charged particle 
accelerators, we express the atomic position in terms of a 
"phase angle" 0, with a periodicity of 2L. Two examples 
of potentials are depicted in Fig. ^ a sinusoidal poten- 
tial and a chained harmonic one. The potential maxima 
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occur at = — 7r/2 + 2mT (n = 0,±1,±2...), while the 
minima are at (j) — 7r/2 + 2mT. Now suppose a particle 
initially located near the potential maximum accelerates 
to the potential minimum for a time r, at which moment 
the original potential is switched off and a new set of 
potential wells identical to the previous ones, but with a 
phase shift of tt, is switched on. Then the particle will 
find itself again near the top of the potential well, and 
continue to be accelerated. 

In order to describe quantitatively this acceleration 
scheme, it is convenient to introduce the concept of the 
synchronous particle. This is a particle that travels ex- 
actly the distance L (or n in terms of phase angle) dur- 
ing the time interval r. We denote its phase and velocity 
right after the switching by and Vs, respectively. By 
definition, then, the initial phase 0sof the synchronous 
particle remains unchanged. 

The kinetic energy gained by the synchronous particle 
during one acceleration stage is given by 

AExiq^s) = V{<j)s) - Vi<j)s + tt). (27) 

We can regard this increase as originating from a contin- 
uously acting, average force 

F(0,) = AEKi(j)s)/L. 

Consider now a non-synchronous particle with initial 
phase (j) = A.(j) + (j)s and velocity v — Av -t- Vs . Its motion 
relative to the synchronous particle is governed approxi- 
mately by the equation of motion 

^A</.--^[i^(</>, + A<^)~F(</>,)]. (28) 

Depending on the initial conditions, Eq. ( p8| ) describes 
either stable or unstable dynamics of the motion of a 
collection of particles. It is easy to see that for stable 
acceleration, we must have 

-7r/2 <(j)s<0. 

Fig. ^ plots the phase stability diagram for the two 
potentials depicted in Fig. ^, obtained by numerically 
solving Eq. (p8|). For a given 0s, the closed curve de- 
picts the boundary of a "bucket" , inside of which a sta- 
ble acceleration is achieved. Since AEk{4is) decreases 
from 2Vb to when (/)s changes from — 7r/2 to 0, we see 
there is always a trade-off between fast acceleration and 
large bucket size, i.e., large acceptance. Furthermore, for 
— 7r/2 < (f)s < — 7r/12, the chained harmonic potential 
gives a larger bucket compared to the sinusoidal poten- 
tial. 

To test our scheme on a condensate, we have solved 
the time-dependent Gross-Pitaevskii equation ( [l0| ) with 
the time-dependent potential of Fig. |[ Our results are 
illustrated in Fig. which shows the velocity of the 
condensate at each of the switching times of the external 
fields, for both harmonic and sinusoidal potentials. The 



switching is chosen to occur at the moment the center 
of mass of the condensate travels a distance L. For this 
calculation, we choose the initial wave function to be a 
Gaussian centered at 00 = — 7r/6 with a width of Q.IL. 
The inserts in the figure show the normalization of the 
wave function, which characterizes the loss of atoms dur- 
ing their acceleration. This is calculated by discarding 
the atoms leaked outside of the well that traps the cen- 
ter of the mass of the condensate. Less loss is observed 
for the harmonic potential than for the sinusoidal one, 
a result of the larger bucket size in the former case (see 
Fig. ^). We also observe a weaker influence of the con- 
densate mean-field energy in the case of a harmonic po- 
tential. However, while atomic collisions do increase the 
losses and decrease acceleration rates, their effect is not 
very dramatic, a result similar to that already encoun- 
tered with the chirped lattice accelerator. In addition, we 
checked the momentum distribution of the wave function 
during the acceleration and found no significant varia- 
tions of the momentum width, hence, monochromaticity 
is also preserved in this scheme. 

The solid line in Fig. |l^ gives the velocity of a classical 
particle under the action of a continuous and constant 
force i^(0o) = AEk{4'o)/L. As the figure shows, the 
average acceleration a of the condensate can indeed be 
approximated by a = F{(j)o)/m. 

We finally note that in the case of alkali atoms, an 
acceleration of ^ 10^ m/s^ requires a well depth of 10 
mK and width L — 100 fim is needed. In view of current 
advances in the technology of magnetic microtraps 
these numbers should not pose as an unsurmountable 
challenge for experimentalists. 

IV. SUMMARY AND CONCLUSION 

In conclusion, we demonstrated two coherent mecha- 
nisms to accelerate Bose-Einstein condensates with very 
little atom loss while preserving its monochromaticity. 
The first mechanism, based on the use of an accelerated 
optical lattice, can be interpreted as resulting from a se- 
quence of Bragg resonances, which can lead to a remark- 
able preservation of the monochromaticity of the atomic 
sample under appropriate conditions. This analysis en- 
abled us to identify an optimal regime of acceleration 
rates for reasonable experimental parameters. The sec- 
ond scheme uses a time-varying potential analogous to 
those used in the synchronous acceleration of charged 
particle. Acceleration is achieved by switching on and 
off the potential in an appropriate way such that the 
condensate keeps gaining kinetic energy at the expense 
of its potential energy. In contrast to the first scheme 
where only quasi-discrete modes are populated, here the 
momentum of the condensate evolves continuously. 

In future work, it will be important to optimize the 
time-dependence of the accelerators, so as to achieve a 
factor of merit as close to unity as possible. We propose 
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to approach this problem using genetic learning algo- 
rithm techniques. It will also be necessary to go past the 
mean-field theory in order to determine the higher-order 
coherence properties of the accelerated beams, in partic- 
ular their atom statistics. In addition, these accelera- 
tion devices are potentially important to generate non- 
classical matter wave fields to be used in improving the 
signal to noise ratio in future atom optical sensors. The 
discrete set of momentum sidemodes that are macroscop- 
ically excited in the sequential Bragg resonances make 
this problem particularly well suited for quantum opti- 
cal "essential modes" approaches familiar from quantum 
nonlinear optics. 



ACKNOWLEDGMENTS 

This work is supported in part by the U.S. Office of 
Naval Research under Contract No. 14-91-J1205, by the 
National Science Foundation under Grant No. PHY98- 
01099, by the U.S. Army Research Office, and by the 
Joint Services Optics Program. 



[1] T. L. Gustavson, P. Bouyer, and M. A. Kasevich, Phys. 

Rev. Lett. 78, 2046 (1997). 
[2] J. Fujita, S. Mitake, and F. Shimizu, Phys. Rev. Lett. 

84, 4027 (2000). 
[3] A. Soroko, J. Phys. B 30, 5621 (1997). 
[4] O. Zobay, E. V. Goldstein, and P. Meystre, Phys. Rev. 

A 60, 3999 (1999) 
[5] D. Jaksch, H.-J. Briegel, J. I. Cirac, C. W. Gardiner, and 

P. Zoller, Phys. Rev. Lett. 82, 1975 (1999). 
[6] G. K. Brennen, C. M. Caves, P. S. Jessen, and I. H. 

Deutsch, Phys. Rev. Lett. 82, 1060 (1999). 
[7] D. Miiller, D. Z. Anderson, R. J. Grow, P. D. D. 

Schwindt, and E. A. Cornell, Phys. Rev. Lett. 83, 5194 

(1999). 

[8] M. Key, I. G. Hughes, W. Rooijakkers, B. E. Sauer, E. A. 

Hinds, D. J. Richardson, and P. G. Kazansky, Phys. Rev. 

Lett. 84, 1371 (2000). 
[9] N. H. Dekker, C. S. Lee, V. Lorent, J. H. Thywissen, S. P. 

Smith, M. Drndic, R. M. Westervelt, and M. Prentiss. 

Phys. Rev. Lett. 84, 1124 (2000). 
[10] W. Hansel, J. Reichel, P. Hommelhoff, and T. W. Hansch, 

Phys. Rev. Lett. 86, 608 (2001). 
[11] E. V. Goldstein and P. Meystre, Phys. Rev. Lett. 80, 

5036 (1998). 

[12] E. V. Goldstein, O. Zobay, and P. Meystre, Phys. Rev. 

A 58, 2373 (1998). 
[13] G. A. Prataviera, E. V. Goldstein, and P. Meystre, Phys. 

Rev. A 60, 4846 (1999). 
[14] P. Meystre, Atom Optics, Springer Verlag, New York 

(2001), in press. 



[15 
[16 

[17 
[18 
[19 
[20 
[21 

[22; 
[23; 

[24; 
[25; 

[26; 

[27 

[28' 
[29; 

[30 
[31 

[32; 

[33 
[34; 
[35; 

[36 

[37 
[38; 
[39 



H. L. Bethlem, G. Berden, A. J. A. van Roij, F. M. H. 
Crompvoets, and G. Meijer, Phys. Rev. Lett. 84, 5744 
(2000). 

J. A. Maddi, T. P. Dinneen, and H. Gould, Phys. Rev. A 
60, 3882 (1999); H. Bethlem, G. Berden, and G. Meijer, 
Phys. Rev. Lett. 83, 1558 (1999); H. Bethlem, G. Berden, 
F. M. H. Crompvoets, R. T. Jongma, A. J. A. van Roij, 
and G. Meijer, Nature 406, 491 (2000). 
M. B. Dahan, E. Peik, J. Reichel, Y. Castin, and C. Sa- 
lomon, Phys. Rev. Lett. 76, 4508 (1996). 
E. Peik, M. Ben Dahan, I. Bouchoule, Y. Castin, and 
C. Salomon, Appl. Phys. B 65, 685 (1997). 
E. Peik, M. B. Dahan, I. Bouchoule, Y. Castin, and C. Sa- 
lomon, Phys. Rev. A 55, 2989 (1997). 
Q. Niu, X.-G. Zhao, G. A. Georgakis, and M. G. Raizen, 
Phys. Rev. Lett. 76, 4504 (1996). 

C. F. Bharucha, K. W. Madison, P. R. Morrow, S. R. 
Wilkinson, B. Sundaram, and M. G. Raizen, Phys. Rev. 
A 55, R857 (1997). 

S. R. Wilkinson, C. F. Bharucha, K. W. Madison, Q. Niu, 
and M. G. Raizen, Phys. Rev. Lett. 76, 4512 (1996). 
K. W. Madison, C. F. Bharucha, P. R. Morrow, S. R. 
Wilkinson, Q. Niu, B. Sundaram, and M. G. Raizen, 
Appl. Phys. B 65, 693 (1997). 

J. Zapata, A. M. Guzman, M. G. Moore, and P. Meystre, 
Phys. Rev. A 63, 023607 (2001). 

W. H. Oskay D. A. Steck, B. G. Klappauf, and M. G. 
Raizen, Laser Phys. 9, 265 (1999). 

D. L. Haycock, P. M. Alsing, I. H. Deutsch, J. Grondalski, 
and P. S. Jessen, Phys. Rev. Lett. 85, 3365 (2000). 

K. Berg-S0rensen and K. M0lmer, Phys. Rev. A 58, 1480 
(1998). 

D. -I. Choi and Q. Niu, Phys. Rev. Lett. 82 2022 (1999). 
B. P. Anderson and M. A. Kasevich, Science 282, 1686 

(1998) . 

J. Stenger, S. Inouye, A. P. Chikkatur, D. M. Stamper- 
Kurn, D. E. Pritchard, and W. Ketterle , Phys. Rev. 
Lett. 82, 4569 (1999). 

E. W. Hagley, L. Deng, M. Kozuma, M. Trippenbach, 
Y. B. Band, M. Edwards, M. Doery, P. S. Julienne, 
K. Helmerson, S. L. Rolston, and W. D. Phillips, Phys. 
Rev. Lett. 83, 3112 (1999). 

L. Deng, E. W. Hagley, J. Wen, M. Trippenbach, Y. B. 
Band, P. S. Julienne, J. E. Simsarian, K. Helmerson, S. L. 
Rolston, and W. D. Phillips, Nature 398, 218 (1999). 
E. V. Goldstein and P. Meystre, Phys. Rev. A 59, 3896 

(1999) . 

M. Trippenbach, Y. B. Band, P. S. Julienne, Phys. Rev. 
A 62, 023608 (2000). 

S. Inouye, T. Pfau, S. Gupta, A. P. Chikkatur, A. Gorlitz, 
D. E. Pritchard, and W. Ketterle, Nature 402 , 641 
(1999). 

M. Kozuma, L. Deng, E. W. Hagley, J. Wen, R. Lutwak, 
K. Helmerson, S. L. Rolston, and W. D. Phillips, Phys. 
Rev. Lett. 82, 871 (1999). 

O. Morsch, J. H. Miiller, M. Cristiani, and E. Arimondo, 



eprint cond-mat/0103466 



P. B. Blakie and R. J. Ballagh, J. Phys. B: At. Mol. Opt. 
Phys. 33, 3961 (2000). 

D. Jaksch, C. Bruder, J. I. Cirac, C. W. Gardiner, and 



7 



p. Zollor, Phys. Rev. Lett. 81, 3108 (1998). 
[40] C. Raman, M. Kohl, R. Onofrio, D. S. Durfee, C. E. 

Kuklcwicz, Z. Hadzibabic, and W. Ketterle, Phys. Rev. 

Lett. 83, 2502 (1999). 
[41] R. Onofrio, C. Raman, J. M. Vogels, J. R. Abo-Shaeer, 

A. P. Chikkatur, and W. Ketterle, Phys. Rev. Lett. 85, 

2228 (2000). 

[42] S. Burger, F. S. Cataliotti, C. Fort, F. Minardi, M. In- 
guscio, M. L. Chiofalo, and M. P. Tosi, to appear in Phys. 
Rev. Lett. 86, 14 May (2001). 

[43] J. Roichol, W. Hansel, P. Hommelhoff, and T. W. Hansch, 
Appl. Phys. B 72, 81 (2001). 




t[ns] 

FIG. 1. Mean displacement of the atomic cloud: The 
simulation uses 10^ Sodium atoms initially at rest with a 
Gaussian spatial distribution of longitudinal width 50 fxm and 
transverse width 5 /nm. The laser wavelength is \l = 985 nm 
and the lattice depth Vq is chosen to bo half the lattice re- 
coil energy. The Sodium mass is m = 3.82 • IQ-^^ kg and 
the s-wave scattering length Os = 4.9 nm. Two upper curves: 
oiat = 724.1 ms~^. Two lower curves: aiat = 327.8 ms~^. 
Solid lines: classical particles, dotted lines: BEC. Unless oth- 
erwise stated, the same parameters are used through Fig. 7. 
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FIG. 2. Mean velocity of the atomic cloud: 

aiat = 724.1 ms~^ (three upper curves), aiat = 327.8 ms^^ 
(three lower curves). Solid lines: classical particles; dotted 
lines: BEC; dashed lines: instantaneous lattice velocity. 
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(a) t=433 




(b) t=433 us 
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FIG. 3. Momentum distribution n(fe) of an ensemble of 10® 
cleissical particles for aiat=(a) 327.8 ms~^, (b) 724.1 ms~^, 
in units of kh- Vertical dotted line: instantaneous lattice 
momentum at time of plot. 
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FIG. 5. Temporal evolution of the mode population for 
an acceleration of 327.8 ms~^. Upper graph: n = to n = 4 
modes. Lower graph: n = — 1 mode. Solid line: solution of 
truncated coupled-mode equations; dotted line: direct solu- 
tion of the Gross-Pitaevskii equation. 
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FIG. 6. Maximum population of the ith mode as a 
function of acceleration, obtained from the truncated cou- 
pled-mode equations. 



FIG. 4. Momentum distribution 0(fc) of a Bose-Einstein 
condensate of W atoms for aiat= (a) 327.8 ms"^ (b) 
724.1 ms~^, in units of fci. Vertical dotted line: instanta- 
neous lattice momentum at time of plot. 
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t [jis] 

FIG. 7. Population of the zcroth momentum mode for 10^ 
(solid), 10^ (dotted) and lO** atoms (dashed) in the conden- 
sate, obtained from a simulation of the full Gross-Pitaevskii 
equation for an acceleration oiat=(a) 327.8 ms~^. 
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FIG. 8. Potential energy as a function of phase angle. 
The dashed line represents a sinusoidal potential, the solid 
line represents a chained harmonic potential. For both cases, 
the well depth is 2Vo. 



FIG. 9. Phase stability diagram for various values of 0s for 
(a) the harmonic potential and (b) the sinusoidal potential. 
The unit for is ,^2nVo / (mL^). 
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FIG. 10. Condensate velocity at the times of the suc- 
cessive switchings for (a) the harmonic potential and (b) 
the sinusoidal potential. Stars: non-interacting condensate 
with {/ = 0; triangles: interacting condensate with U = 10. 
For practical choices of parameters, the typical value for U 
in our dimensionless units is on the order of 1 to 10. In- 
serts: normalization of the wavepackcts. Time in units of 
2mL^ /{hn^), velocity in units of h-K/{mL), and energies in 
units of h^Tv^ /{2mL^). In this example, Vb = 150. (In prac- 
tice, Vb should have a much larger value for realistic choices 
of parameters. However, in our calculations large values of Vb 
gives rise to too much numerical noise.) 
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